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In this paper we present the constraints on cold dark matter (CDM) isocurvature contributions 
to the cosmological perturbations. By employing Markov Chain Monte Carlo method (MCMC), 
we perform a global analysis for cosmological parameters using the latest astronomical data, such 
as 7-year Wilkinson Microwave Anisotropy Probe (WMAP7) observations, matter power spectrum 
from the Sloan Digital Sky Survey (SDSS) luminous red galaxies (LRG), and "Union2" type la 
Supernovae (SNIa) sample. We find that the correlated mixture of adiabatic and isocurvature 
modes are mildly better fitting to the current data than the pure adiabatic ones, with the minimal 
given by the likelihood analysis being reduced by 3.5. We also obtain a tight limit on the fraction 
of the CDM isocurvature contributions, which should be less than 14.6% at 95% confidence level. 
With the presence of the isocurvature modes, the adiabatic spectral index becomes slightly bigger, 
nf^^ — 0.972 ±0.014 (1 a), and the tilt for isocurvature spectrum could be large, namely, the best fit 
value is nf° = 3.020. Finally, we discuss the effect on WMAP normalization priors, shift parameter 
R, acoustic scale I a and z*, from the CDM isocurvaure perturbation. By fitting the mixed initial 
condition to the combined data, we find the mean values of R, Ia and z, can be changed about 
2.9cr, 2.8a and 1.5cr respectively, comparing with those obtained in the pure adiabatic condition. 



I. INTRODUCTION 

The accumulation of WMAP seven year measurement on cosmic microwave background radiation (CMB) [l|, 
associated with observations of SDSS , provide wealthy information on the anisotropics and inhomogeneities of our 
universe, in light of which the perturbation theory has been tested in certain level. Currently, various observational 
data favor the simplest concordance cosmological model, which has six free parameters and the pure adiabatic initial 
condition [l|, [^-Q . Although the concordance model fits the data quite well, it is always worthy to consider alternative 
candidates. Namely, it is important to study observational constraints on initial states of cosmological perturbations 
at reheating surface. 

Generically, there exists two classes of modes of cosmological perturbations, with one being adiabatic of which the 
trajectory is parallel to the background evolution, while the other isocurvature of which the trajectory is orthogonal to 
the background evolution [7l-[T5| . However, a first lesson from observational data is that the primordial fluctuations are 
nearly adiabatic; additionally, an isocurvature mode is also expected to be negligible as is predicted by the simplest 
inflation model in terms of a single infiaton field and a rapidly reheating process [l^. Therefore, one usually makes 
the data fitting without considering the isocurvature mode. 

Accompanied with developments of inflationary cosmology, models of multiple field inflation were extensively stud- 
ied in the literature p7l - [27j . which predicted an existence of primordial isocurvature fluctuations. These primordial 
isocurvature modes could be transferred into cosmological perturbations after reheating, such as Bayon, CDM, DE, 
and neutrino respectively [28l - l3]| . Consequently, the attention on isocurvature modes has been awaken in recent 
years miH. 

One may notice that, although the pure isocurvature primordial perturbation has been ruled out by the Boomerang 
and MAXIMA- 1 data [sj] already, a mixture of adiabatic and isocurvature modes can be in agreement with the 
current data fortunately. In the literature [35l - l46j . the studies on constraining the isocurvature fluctuation have 
been performed extensively using various observational data, such as CMB, LSS, integrated Sachs- Wolfe effect or 
Lyman-a forest data. With the WMAP7 data, the totally un-correlated and anti-correlatcd adiabatic, non-adiabatic 
perturbations are constrained which arc performed by flxing the correlation coefficients to be or —1, respectively. 
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In this paper, we study the constraints on the isocurvature modes of cosmological perturbations in hght of the 
latest observational data. We consider a generic scenario that the adiabatic modes and isocurvature ones are allowed 
mixed through a correlation matrix which could be arising from a time- varying background trajectory. In the detailed 
analysis, we treat the coefficients of correlation matrix as free parameter and make the data fitting. Comparing with 
the previous results in the literature, our result shows a slight improvement on the final constraints of the mixed 
initial condition parameters. Namely, CDM isocurvature components are stringently limited, the contribution from 
the isocurvature modes arc only allowed in small scales, also the error bars of initial condition parameters become 
smaller than the past works in the literature, which can be observed from the one dimensional probability distribution 
of the fraction parameter a as analyzed in the main context. 

The WAMP normalization priors, R, Ia and encoding the information of background cosmic distances, can 
be applied to greatly simplify the numerical calculations of determining cosmological parameters, such as the EoS 
of dark energy. It was found that these priors could be sensitive to the peak locations and local structures of the 
CMB temperature power spectrum [43|. However, it is well-known that these quantities can be affected by CDM 
isocurvature perturbation. Therefore, we study the effects on the WMAP normalization priors given by the WMAP 
group from the isocurvature mode perturbation. 

The outline of this paper is as follows. In Section II, we describe the parametrization of cosmological perturbations 
with adiabatic and isocurvature modes mixed, and then we consider a specific example to illustrate the effects of 
isocurvature perturbation imprinted on the CMB temperature power spectrum and LSS matter power spectrum 
respectively. In Section III, we perform a global analysis and introduce the data we applied and the parameters used 
in the data fitting. The constraints on these parameters are present in detail in Section IV, and the corresponding 
effects on reduced distance parameters arc discussed. Finally, Section V includes the conclusions and discussion. 



II. PARAMETRIZATION OF CORRELATED ADIABATIC AND ISOCURVATURE PERTURBATIONS 

For model-independent consideration, the initial conditions for the correlated mixture of adiabatic and isocurvature 
modes can be written as 

k 



where fcg is the pivot scale. Here, A^^ and n^J are 2x2 symmetric matrices which characterize the amplitude and 
power spectrum index, respectively. We have 



An = f Af' v/If^4r cosA^, 

' cos A A'^° ' ' ^ ' 



where cos A = ^adi.iso^ ^ ^adi^iso describes the correlation between the adiabatic mode and the isocurvature mode 
[20| . ^f'^' and are the amplitude of adiabatic and isocurvature modes, respectively. The index matrix ri^J is given 
by: 



^cor ^iso 



(3) 



where ng"^' and are the spectrum indices for adiabatic and isocurvaure modes. Following Ref. [s^, we assume 
n™'- = (nf^ + n'/°)/2 for simphcity. 

Both adiabatic and isocurvature perturbations seed the structure formation of our universe, and can be imprinted in 
CMB temperature, polarization power spectrum, as well as matter power spectrum of galaxies surveys. Symbolically, 
the anisotropics of CMB photons can be written as 



6T f 6T\ f 6T\ f 5T 



(4) 



Consequently, the CMB power spectrum is given by 



C/ ' = Af'Cf + A^Cr + 2^ i^isocf cos A , (5) 
where initial power spectrum can be obtained as follows. 



3 



with i, j standing for adiabatic or isocurvaturc and 0^ are the transfer functions of the radiation at initial mo- 
ment. There are similar formulas for CMB polarization power spectra Cf^, Cf^ and temperature-polarization cross 
correlation power spectrum Cj^. 

Moreover, the matter power spectrum P{k) can be written as: 

P{k) = Af'P'"^\k) + A'/°P'^°(/fc) + 2^JAf^Af°P'''^'^''°{k) cos A , (7) 

with 

PHk)^i^-j T\k)T\k), (8) 
where T^{k) are the transfer functions of matter component at initial moment. 





FIG. 1: The CMB temperature power spectra and LSS matter power spectra are plotted using different initial conditions in 
the upper and lower panels, respectively. The red solid lines are given by the best fit model with the mixed initial condition. 
The blue dash lines and the black dash-dot lines are given by choosing = 1 x 10"^" and cos A = —1, respectively. 

In order to show the effect of isocurvature perturbations, in FiglT] we plot the CMB temperature power spectra 
(upper panels) and LSS matter power spectra (lower panel). The red solid lines are given by the best fit model with 
mixed initial conditions by fitting the current observational data: = 0.0226, uJc = 0.1144, A'f' = 2.407 x 10"^, 
Ajf° = 3.885 X 10-1^ nf^ = 0.968, n'^° = 3.020, cos A = 0.149. The blue dash line and black dash-dot hues are 
given by choosing Af° = 1 x 10^^" and cos A = — 1 instead, respectively. One can see that by taking into account 
the CDM isocurvaure modes and its correlation with adiabatic perturbations, the hight and location of the peaks of 
TT spectra will be modified, and the amplitude of matter power spectrum on small scales arc changed significantly. 
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III. OBSERVATIONAL CONSTRAINTS 



A. Data 



We perform a global analysis by employing a modified MCMC package CosmoMC ^ [i^ , in which we have made an 
extension of the initial perturbations by including correlated adiabatic and CDM isocurvature modes. In computation 
of CMB we have included the WMAP7 temperature and polarization power spectra with the routine for computing 
the likelihood supplied by the WMAP team^ Since the CDM isocurvature perturbations affect the power spectrum 
on small scales, we also include some small-scale temperature anisotropics measurements, ACBAR [i^, CBI [H^l 
and Boomerang [sH - Furthermore, we use the matter power spectrum measured by the observations of luminous 
red galaxie (LRG) from SDSS data release seven 0, and the "Union 11" supcrnovae datasets calibrated by SALT2 
template [521 ■ 



B. Parameters 



As pointed out in Refs.[39|, |4^, the constraints of the spectral index may depend on the pivot scale. In order to 
avoid such dependence, we do not take the spectra indices n^^ as free parameters in our analysis. We take two different 
amplitudes As{ki) and As{k2) at two pivot scales ki = 0.002 and k2 = 0.05 to be free. And the spectra indices can 
be derived by 

_ ^ ^ in[pHk2)/rHk,)] ^ in[Ay(fc2)My(fcO] 

ln[fc2/fci] ln[fc2/fci] 
We assume that the overall amplitude is composed of two components: 

A, = Af' + A^^'° , (10) 

and the fraction of adiabatic part is defined by 

a = Af'/A,. (11) 

Thus, our most general parameter space is 

P = {ujb,ujc,Qs,T,cosA,ai,a2,As{ki),As{k2)}. (12) 

where uJb = ^bh^, = Of,/i^ denotes physical baryon density and cold dark matter density relative to critical density 
respectively; Qg = lOOfj is the ratio of sound horizon to angular diameter distance at decoupling, t characterizes the 
optical depth to reionization. ^^(fci), ^3(^2), ai and 012 are the amplitudes and the fractions of adiabatic mode at 
two scales, respectively. In our calculation, we assume the flat universe and dark energy is the cosmological constant. 



TABLE I: The mean values and 1 a errors of As given by EQ.fTUlat two different pivot scale, 02 and cos A, while the limit of 
ai is for 95% confidence level. 



log[10"M4fci)] 


log[10^'M,(fc2)] 


Otl 


a2 


COS A 
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TABLE II: The mean values and 1 a errors of the parameters related to the mixed initial adiabatic and CDM isocurvature 
perturbations. 





log[10"Mf'] 


nf 
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3.204 ± 0.036 
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3.165 ± 0.044 


0.972 ±0.014 
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FIG. 2: One dimensional probability distributions of initial conditions associated parameters are plotted with red solid lines. 
For logflO^'' AJ"*'] and nj''', we specially make a comparison with the results obtained from the pure adiabatic initial conditions 
which are plotted with black dash lines. We also plot the adiabatic perturbation modes fraction parameter oi and 02 using 
the purple dash-dot line and blue solid line, respectively. 

IV. GLOBAL FITTING RESULTS 

In this section, we present the main constraints on the cosmological parameters, as well as the parameters related 
to the CDM isocurvature perturbations. Then, we discuss the effects on the WMAP normalization priors when taking 
into account CDM isocurvature perturbations. 

A. Constraints on cosmological parameters 

In table I and tablcjll] we list the constraints on the initial condition associated parameters by fitting the latest 
observational data. The parameters presented in table I are the free parameters during our calculation, while those 
in table HI] are the derived parameters from the those in table I according to EQ. (| WTT|) . and both results are given 
by the statistics of the samples given by our MCMC approach. From the mixed initial conditions, the constraints on 
the adiabatic perturbation are \og[10^° Af'] = 3.165 ± 0.044 and nf' = 0.972 ± 0.014 at 68% confidence level. The 
adiabatic spectrum index becomes larger, when comparing with the constraint nf^^ = 0.964 ± 0.010 (Icr) from the 
pure adiabatic initial condition. For the CDM isocurvature mode, the constraints are log[10^°A!f°] = ~1.375'ti'^oe 
and nf° = 2.246toAlt (68% C.L.). Furthermore, the coefficient of the correlation is limited as -0.113 < cos A < 0.473 
at 95% C.L.. Comparing with the minimal of these two cases, we find that a weakly positive-correlated mixture of 
CDM isocurvature and adiabatic perturbations is mildly favored by the current data. However, the fraction of CDM 
isocurvature mode can not be larger than 14.9% at 95% confidence level. The constraining power mainly comes from 
WMAP7 temperature and polarization power spectra, since too large fraction of the CDM isocurvature mode will 
lead to the overstated modification of the shape of acoustic peaks in TT and TE power spectra, as illustrated in FiglT] 
Our results are consistent with the previous results reported in Ref . [31 • 

In Figl^l we plot one dimensional posterior probability distributions of A^'^', Af°, nf^' and nf°, as well as cos A 
and a at the two different pivot scales. In Fig|3]we also plot the one dimensional probability distributions of the other 
cosmological parameters obtained from the mixed initial conditions and pure adiabatic condition cases. The red solid 
lines are given by the mixed initial conditions, while the black dash lines are from the pure adiabatic one. One can 
see that with the mixed initial conditions, the results favor a slightly smaller and larger 9 due to the correlation of 
the CDM isocurvature perturbations with other cosmological parameters when considering the non-adiabatic initial 
condition in our analysis. 

We find that there exist apparent correlations among the parameters related to initial conditions and other back- 
ground cosmological parameters, as illustrated in Fig. 21 The coefficient cos A is correlated with fi™, since a positive 
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FIG. 3: One dimensional probability distributions of some background cosmological parameters, the red solid lines are given 
by the mixed initial condition, while the black dash-dot lines are from the pure adiabatic conditions. 




FIG. 4: Two dimensional contours among some cosmological parameters. 
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correlated adiabatic and isocurvature component raise the peaks of TT spectrum which can be compensated by lower 
iljn- A strongly negative correlation between nf° and Af° is shown in the right upper panel of Fig. 2] 



B. Information on reduced distance parameters 



TABLE III: Median 1 a constraints on WMAP normalization priors using full WMAP7 data for different initial perturbations. 





R 


Ia 




2» 


WMAP only 


Adiabatic 


-1 79-|+U.Ui8 
^- '^^-0.019 


302.3 ±0.7 


146.7 ± 1.5 


1091.1 ±0.9 


Mixed 


1. (U1_Q Q28 


301.5+i-; 


147.3 ± 1.9 


1089.8 ± 1.4 


Full data 


Adiabatic 


1.729 iO.OfO 


301.8 ±0.6 


145.7 ±0.9 


1091.3 ±0.6 


Mixed 


i. /UU_Q (J25 


300.1 ± 1.3 


146.4 ± 1.5 


1090.4 ± 1.0 



The WMAP normalization priors given by WMAP group include the "shift parameter" R, the "acoustic scale" Ia 
and the photon decoupling epoch . R and I a correspond to the ratio of angular diameter distance to the decoupling 
era over Hubble horizon and sound horizon at decoupling respectively, given by 



R{z,) = ^n„,H^r{z,) , (13) 

Ia{z*) = 7rx(z*)/rs(2:*) , (14) 

where r(z*) and 7's(z*) denote the comoving distance to and comoving sound horizon at respectively. The 
decoupling epoch is given by 53 1 



z* = 1048[1 ± 0.00124:{nbh^)-° ™][l + .gi(an/i^)^'] , (15) 

where 

_ 0.0783(»b/i^)-°-^38 0.560 

" l±39.5(r!b/i2)0-763' 92 - i + 21.1(r!6/i2)i-8i- ^ ^' 

The WMAP normalization priors encode in part of the CMB information and can be used to constrain cosmo- 
logical parameters to some extent. They are derived parameters from the CMB power spectra based on the fiducial 
cosmological model. Thus, they arc model dependent [47l. l54i|. By using the WMAP7 data only, in Fig. [5] we show 
the one-dimensional probability distributions of the WAMP normalization priors in two different initial perturbation 
conditions. The red solid lines are given by pure adiabatic initial perturbations, while the black dash-dot lines are 
with the mixed case. Wc find the obvious difference in the probability distributions of R, I a, and z*. As we know, 
the WMAP normalization priors mainly include the information on the oscillatory structures of the CMB power 
spectrum. As shown in FiglTJ the TT power spectrum can be modified by the CDM isocurvature perturbations, 
namely the locations and the height of the peaks are modified when taking into account the contributions of CDM 
isocurvature perturbations. And these effects can lead to the difference of the mean values and errors of WMAP 
normalization priors. 

In table Hill wc list la constraints on WMAP normalization priors from the full WMAP7 data for different initial 
conditions. One can sec that, by taking into account the contribution from the initial CDM isocurvature modes, 
constraints of the WMAP normalization priors are obviously changed. The mean values of R is changed about la 
when comparing with the pure adiabatic case, while for I a and z^, the changes are about la and lAa, respectively. 

Furthermore, we also provide the 1 a constraints on R, I a, z* and using the combined data sample of WMAP7, 
ACBAR, CBI, Boomerang, as well as SDSS LRG and SN la in table Hill As we expect, combining the data of SN 
and LSS can help in constraining the other background cosmological parameters and shrink error bars of R, I a, -z* 
and significantly. When considering the CDM isocurvature mode perturbation, we find the differences between the 
constraints from pure adiabatic perturbation and mixture of adiabatic and isocurvature perturbations are improved 
obviously. In this case, the difference of the mean values of R, Ia, z* and become 2.9a, 2.8a, 1.5a and 0.7cr, 
respectively. 
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FIG. 5: One dimensional posterior distributions of I a, R and with WMAP7 data, the red solid lines are given by the mixed 
initial perturbations, while the black dash lines are from pure adiabatic perturbations. 



V. SUMMARY 



In this paper, we study the constraints on mixed adiabatic and isocurvature modes of cosmological perturbations. 
Using the current observational data, such as WMAP7 CMB power spectrum, matter power spectrum of SDSS data 
released seven LRG data and SNIa "union2" sample, we find an adiabatic initial condition with the presence of certain 
isocurvature modes can explain the experiments better than a pure adiabatic initial condition. Moreover, we obtained 
more stringent constraints on the parameters of cosmological perturbations by virtue of the improvement of accuracy 
of observational data in recent years. Our result shows that the spectral index of isocurvature perturbation has a 
very blue tilt. This quantity may lead to either enhancement or depression of power spectrum on small scales, which 
depends on the correlation angle of isocurvature and adiabatic modes. 

Given the WMAP normalization priors are widely used in performing testing of cosmological models, we provide 
the comparison on constraints of i?, Ia and from different initial conditions. Since WMAP normalization priors 
are derived parameters from the CMB power spectra, the constraints are changed obviously when the isocurvature 
mode modifies of the shape of the peaks and troughs in CMB power spectra. 

As an end, we would like to mention that, from the point of view of information criteria, detailed constraints 
on cosmological parameters depend on the set of parameters of the model which is compared with observational 
datajssj. The simplest model with purely adiabatic scale-invariant primordial power spectrum is able to capture 
the most relevant clues of early universe physics, however, more information deserves to be explored along with the 
improvement of observational data from forthcoming experiments, we expect the global analysis on both CMB power 
spectrum and matter power spectrum on small scales will be crucial to constrain the isocurvature perturbation. 
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